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Diluting the inflationary axion fluctuation by a stronger QCD in the early Universe 
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We propose a new mechanism to suppress the axion isocurvature perturbation, while producing 
the right amount of axion dark matter, within the framework of supersymmetric axion models 
with the axion scale induced by supersymmetry breaking. The mechanism involves an intermediate 
phase transition to generate the Higgs /x-parameter, before which the weak scale is comparable to 
the axion scale and the resulting stronger QCD yields an axion mass heavier than the Hubble scale 
over a certain period. Combined with that the Hubble-induced axion scale during the primordial 
inflation is well above the intermediate axion scale at present, the stronger QCD in the early Universe 
suppresses the axion fluctuation to be small enough even when the inflationary Hubble scale saturates 
the current upper bound, while generating an axion misalignment angle of order unity. 


The non-observation of the neutron EDM requires the 
CP violating QCD angle to be as tiny as \9\ < 10” 10 , 
causing the strong CP problem. An appealing solution of 
this puzzle is to introduce a spontaneously broken global 
Peccei-Quinn (PQ) symmetry [l|. Then 0 corresponds 
to the vacuum value of the associated Nambu-Goldstone 
boson, the axion, which is determined to be vanishing by 
the low energy QCD dynamics Hj. 

An interesting consequence of this solution is that ax- 
ions can explain the dark matter in our universe. Yet, 
the prospect for axion dark matter depends on the cos¬ 
mological history of the PQ phase transition. A possible 
scenario is that the spontaneous PQ breaking occurs af¬ 
ter the primordial inflation is over. In such a case, the 
model is constrained to have the domain-wall number 
-Ydw = 1 , where IVdw corresponds to the integer-valued 
U(1)pqxSU(3) c xSU(3) c anomaly coefficient. Then ax- 
ions are produced mainly by the annihilations of ax- 
ionic strings and domain-walls, which would result in the 
right amount of axion dark matter for the axion scale 
/ a ~ 5 x 10 10 GeV Q. However it appears to be difficult 
to realize this scenario within the framework of a funda¬ 
mental theory such as string theory, since it requires a 
PQ symmetry with V DW = 1, as well as a restored PQ 
phase until some moment after the primordial inflation. 

Another scenario which we will focus on in this paper is 
that U(1)pq is spontaneously broken during the primor¬ 
dial inflation and never restored afterwards. Then the 
model is not subject to the condition Np ,w = 1 , but is 
constrained by the axion isocurvature perturbation ■ 
For instance, from the observed CMB power spectrum, 
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where 9 m ; s and 59 denote the average misalignment angle 
and the angle fluctuation, respectively, for the axion field 
right before the conventional QCD phase transition when 
m a {tQCT>) ~ H(tQcn) with a temperature T(£q C d) ~ 
1 GeV. The relic axion density is given by 
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with Qdm ~ 0.24 being the total dark matter fraction. 
Here we have assumed that \80\ -C |0 m is| an d there is 
no significant evolution of f a from £qcd to the present 
time to so that f a (tQcu) ~ faifo)- In inflationary cos¬ 
mology, the primordial quantum fluctuation of the axion 
field results in 


69 = <56»(£qcd) = 7 69 ( h ) = 7 ^ (3) 

where f a (ti) and H(ti) denote the axion scale and the 
Hubble parameter, respectively, during the primordial in¬ 
flation epoch £/, and the factor 7 is introduced to take 
into account the evolution of 59 from £/ to £qcd- Note 
that the inflationary Hubble scale H(tj) is bounded by 
the tensor-to-scalar ratio of the CMB perturbation as 


r ~ 0.16 
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and the weak gravity conjecture [sj suggests that generic 
axion scales are bounded as 
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where Mpi ~ 2.4 x 10 18 GeV is the reduced Planck mass. 
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To discuss the implication of the isocurvature con¬ 
straint o, one needs to specify the cosmological evo¬ 
lution of the axion scale after the primordial inflation is 
over. If f a (tj ) ~ / a (t 0 ) as has been assumed in most 
of the previous studies, it requires that either H(ti) is 
smaller than its upper bound ~ 10 14 GeV by at least five 
orders of magnitude, so that the CMB tensor mode is too 
small to be observable, or 59 should experience a large 
suppression after the primordial inflation, which appears 
to be difficult to be implemented. 

The above observation suggests a more attractive sce¬ 
nario realizing f a (ti) > /a(to) i in a natural man¬ 
ner. Indeed supersymmetric axion models offer a nat¬ 
ural scheme to realize such a scenario, generating the 
axion scale through the competition between the tachy- 
onic SUSY breaking mass term and a supersymmetric, 
but Planck-scale-suppressed higher dimensional term in 
the scalar potential 10|~13|. One then finds 
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fa (t/) 
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which explains elegantly the origin of an intermediate 
axion scale at present, while giving a Hubble-induced in¬ 
flationary axion scale well above the present axion scale, 
if the supersymmetry (SUSY) breaking mass msuSY at 
present is around TeV scale. Furthermore, this type of 
axion models can be successfully embedded into string 
theory. Specifically, they can be identified as a low en¬ 
ergy limit of string models involving an anomalous U(l ),4 
gauge symmetry with vanishing Fayet-Illiopoulos term 


IZllM- I 11 such string models, the U(l ),4 gauge boson 


is decoupled from the low energy world by receiving a 
heavy mass Ma ~ g 2 Mpi/8-K 2 through the Stiickelberg 
mechanism, while leaving the global part of U(1 )a as 
an unbroken PQ symmetry in the supersymmetric limit. 
Once SUSY breaking is introduced properly, in both the 
present Universe and the inflationary early Universe, the 
residual PQ symmetry can be spontaneously broken to 
generate the axion scales as 0. 


In this paper, we discuss a novel mechanism to sup¬ 
press the axion isocurvature perturbation, while produc¬ 
ing the right amount of axion dark mater, within the 
framework of supersymmetric axion models with the ax¬ 
ion scales given by d 6 j). The isocurvature constraint m 
and the relic axion density 0 suggest that for H(ti) 
near the current upper bound ~ 10 14 GeV, the al¬ 
lowed amount of axion dark matter is maximal when 
f a (to) ~ 10 11 -10 13 GeV, while f a (ti) nearly saturates 
the weak gravity bound (0, e.g. f a (fi) ~ 10 16 - 10 17 GeV. 
Interestingly, the axions scales generated by SUSY break¬ 
ing as © automatically realize such pattern if msusY is 
around TeV scale. More specifically, for the case 


(7) 


the isocurvature bound m reads off 

H(tj) /a08\ 2 /Udm\ ( f a (t 0 ) \ 08 / ITeV \ 

10 14 GeV V 7 / Ua j \ 10 12 GeV/ \m S vsY) ’ 

when combined with ©• This implies that a high 
scale inflation scenario with H(ti ) ~ 10 13 10 14 GeV, 
which would give an observable tensor-to-scalar ratio 
r = 0(0.1-0.01) in the CMB perturbation, can be com¬ 
patible with the axion dark matter Q a = Udm, if the 
axion field fluctuation experiences just a mild suppres¬ 
sion after Q, e.g. 7 = 0 ( 0 . 1 - 0 . 01 ) in 0 . 

To suppress 59 through its cosmological evolution, one 
needs a period with m a (i) > H (t) well before £qcd • On 
the other hand, usually this is not easy to be realized 
because the axion mass should be generated mostly by 
the QCD anomaly in order for the stron g C P problem 
solved by the PQ mechanism. (See Refs. Il5l - ll8j for an 
alternative possibility.) In the following, we propose a 
simple scheme to achieve such a cosmological period by 
having a phase of stronger QCD in the early Universe. 

Our scheme is based on a phase transition at t = 
t/, which will be called the /z-transition in the following 
as it generates the Higgs /i-parameter through the super¬ 
potential term [l9l |. 

li(X)H u H d = KlX ^ Hd , (8) 

where X is a PQ-charged gauge-singlet superfield. 
Specifically, 

X(t<t tl ) = 0 , X(t > tp ) ~ \ J tosusy Mpi , (9) 

so that 

£t(i < tp) = 0 , fi(t > ty) ~ tosusy- ( 10 ) 

With this transition, the weak scale experiences an un¬ 
usual evolution in a way that the weak scale before the 
/i-transition is comparable to the axion scale 0 , as will 
be discussed below. 

To proceed, let us discuss first the key features of the 
scheme, and later present an explicit model to realize the 
whole ingredients. Including the Hubble-induced contri¬ 
bution, the mass of the D-flat Higgs direction H u Hd is 
generically given by 

ml = c^H 2 + ^mlusY + 2|/z | 2 (</> 2 = H u H d ), (11) 

where and are model-dependent parameters of 
order unity. In our scheme, both c$ and are as¬ 
sumed to be negative, so to^, < 0 before the /z-transition. 
Then <f = \JH u H d is stabilized by the competition be¬ 
tween the tachyonic ml\c/)\ 2 and a supersymmetric term 
of 0 (|</>| 6 /Mj,,) in the scalar potential, which results in 

fa(tl) ~ ~ y/H(ti)Mpi, 

fa(ty) ~ ~ \J tosusy Afp/- 


fa(to)/fa{tl) w \JmsusY/H(ti), 


( 12 ) 


















3 


On the other hand, after the /t-transition, > 0 due to 
y ~ tosusy- The resulting weak scale and axion scale at 
present are given by 

<j)(t Q ) = 0(100) GeV, 

fa(t 0 ) ~ X(t 0 ) ~ yjmsusyMpi. (13) 


A simple consequence of the above evolution of H u H<i 
is that the weak scale is comparable to the axion scale 
before the /x-transition: 

4>= tfj) ~ f a = faif < tff). (14) 


This results in a higher QCD scale, i.e. a stronger QCD, 
and therefore a heavier axion mass which might be even 
bigger than the Hubble scale for a certain period. Let 
us estimate the QCD scale Aqcd before the /x-transition, 
which is defined as the scale where the 1-loop QCD cou¬ 
pling blows up, as well as the resulting axion mass rh a . 
For the case with Aqcd < "iTigijhg) < 10 ~ 5 4>, where ihg 
denotes the gluino mass before the /x-transition, we find 


A, 


QCD 


23TeV 


2/11 


30 TeV ) 

tan/A 3/,n / cj) \ 

10 ) ^lO 12 GeV j 


6/11 


(15) 


where tan/3 = ( H u )/(H<i ) at present, and rhg/g 2 (fhg) ~ 
rn g/9z( m g) f° r the gluino mass rrig at present. Here we 
assume that gl(M GVT ) = gl{M GVT ) and y q (M GVT ) = 
y q {M G ut) for the QCD coupling and the quark Yukawa 
couplings. When the temperature T < Aqcd> the axion 
mass during the period of stronger QCD is estimated to 
be 


m a « AQ CD //a- (16) 

On the other hand, if fhg{ A np ) < Aqcd < 10 _5 </>, the 
resulting QCD scale is estimated as 

Aqcd „ -Tev(^) 1 / '( i 5 4 v) W .( 17 ) 

with the axion mass 

rh a « ml /2 A^c D /fa- (18) 


conventional QCD phase transition. As is well known, 
in such a period the axion field experiences a damped 
oscillation, with an amplitude a (averaged over each os¬ 
cillation period) evolving as 

a oc i?“ 3/2 (t), (19) 


where R(t) is the scale factor of the expanding universe. 
Then the spatially averaged vacuum value of the axion 
field is settled down at the minimum of the axion poten¬ 
tial induced by the stronger QCD, while the axion angle 
fluctuation is diluted according to 
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27 Tfaiti) 


T(tg) \ 3/2 H(t T ) 
T(ti) ) 27r/ 0 (£/) ’ 


( 20 ) 


where t = ti denotes the moment when the damped ax¬ 
ion oscillation begins, and t = is the moment when 
it is over. Note that, after the /x-transition, the weak 
scale and the QCD scale quickly roll down to the present 
values, so the axion mass becomes negligible compared 
to H(t) until t ~ £qcd when the Universe undergoes the 
conventional QCD phase transition. Also, the minimum 
of the axion potential induced by the stronger QCD is 
generically different from the minimum of the axion po¬ 
tential at present. As a result, our scheme generates an 
axion misalignment angle of order unity: 



together with an intermediate axion scale at present, so 
gives rise to Q a = Dqm in a natural way. 

In our case, the damped axion oscillation induced by 
the stronger QCD begins at a temperature T(tQ ~ Aqcd 
as fh a is highly suppressed by thermal effects for T 
Aqcd- On the other hand, the scalar field X generat¬ 
ing y through (J5]) is trapped at the origin by thermal 
effects until the Universe cools down to a temperature 
T{tff) ~ tosusy- In fact, our scheme involves a variety 
of dimensionless parameters which affect the naive esti¬ 
mate of the involved scales. We find that there is a large 
fraction of the natural parameter region where the axion 
mass 


m a 


0.4 MeV 


fa \ 

10 12 GeV J 


I Aqcd 
y20TeV 


( 22 ) 


is larger than the Hubble scale 


Here A np denotes the scale where the stronger QCD be¬ 
comes nonperturbative, i.e. around = 8tt 2 /N c with 
N c = 3. Note that the axion potential for the axion mass 
(fT8|) can be obtained by a single insertion of the SUSY 
breaking spurion rhg6 2 to the nonperturbative superpo¬ 
tential Wnp ~ Aqcd induced by the gluino condensation. 

If the stronger QCD scale Aqcd is high enough, there 
could be a period with fh a (t ) > H(t) well before the 


Hfa)c0.2Mev( 1TeVx ^; i2GeV ), (23) 

over the period U < t < with a temperature ratio: 

T{tg)/T(ti) = ^(lO-^lO- 2 ). (24) 

Then the resulting 68 given by (EU1) can be small enough 
to satisfy the isocurvature bound © even when H(tj ) 
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FIG. 1: Upper bound on the inflationary Hubble scale con¬ 
sistent with the axion dark matter, Q a = Udm. Here we 
have taken trig = 3 TeV, tan /3 = 10, and T(Q) = 1 TeV. 
The shaded region is excluded by the Planck results. The 
black solid line is the constraint in the conventional scenario 
with f a (ti) = fa (to)- The magenta lines are for the scenario 
with f a (ti)/f a (t 0 ) = i/H r (t/)/msuSY, but without a stronger 
QCD. The blue lines are for our scheme which leads to a 
further suppression of 56 by the stronger QCD. The SUSY 
breaking mass has been taken msuSY = 1 TeV for the solid 
lines and 10 TeV for the dotted lines. 


saturates its upper bound ~ 10 14 GeV. Note that during 
U<t< tp, 

<f>(t) ~ <P(to) ~ X(t) — X(to) ~ VmsuSYMpi, 

so the corresponding vacuum energy density Vo = 
0( m S\JSY fa(to))- This means that in this period the 
Universe is dominated by the vacuum energy density with 
the Hubble scale give n by (|23|) . which is often called the 
thermal inflation [201 ]. 

It should be stressed that in our scheme the ax¬ 
ion isocurvature perturbation is suppressed by two 
steps. The first suppression is due to f a {to)/faifi) ~ 
\/ m SUSY /H(ti) -C 1, and the second is due to the 
stronger QCD dynamics before the /i-transition, yield¬ 
ing a further suppression by 7 ~ (tosusy/Aqcd) 3 ^ 2 - To 
illustrate the result, we depict in Fig. Q] the upper bound 
on the inflationary Hubble scale H(ti) resulting from the 
isocurvature constraint m for Q a = Qdm- To make a 
comparison, we depict the results for three distinct cases: 
(*) the conventional scenario of f a (ti) = f a (to ) with¬ 
out a stronger QCD, (ii) a scheme with f a (ti)/f a (to ) ~ 
yj H (ti ) /tosusy , but without a stronger QCD, (in) our 
scheme with f a (t 1 ) / f a (to) ~ 'jH(t I )/m SU sy and a 
stronger QCD before the ^.-transition. 

Let us now present an explicit model implementing the 
mechanisms discussed above. As a simple example, we 


consider a model with the following superpotential, 

W = (MSSM Yukawa terms) + AYd "!* 0 
Kl X 2 H u H d , k 2 XY 3 , ac 3 (H u H d )(LH u ) 

M Pl M Pl M P i ' { ] 

where X and Y are PQ-charged gauge singlets responsi¬ 
ble for the /i-transition, L is the MSSM lepton doublet, 
and $ + < 1 > C are I/(l)y-charged exotic matter fields in¬ 
troduced to give a thermal mass to Y. Then the scalar 
potential for the /i-transition is given by 

Vi = m 2 x \X \ 2 + m 2 Y \Y \ 2 + (^-XY 3 +h.c)j 

+ XW (m 6 + 9m 2 |Y| 4 ), (26) 

1V1 PI 

where 

m 2 x = c x H 2 + £xTOg USY + 4|/iA'| 2 ; 

my = cyH 2 + (ymg UgY + otyT 2 , (27) 

for fix = «i H u H d /M P i. Here c X} yH 2 are the Hubble- 
induced masses, CyyTO 2 usy are the SUSY breaking 
masses at zero temperature, and ayT 2 is the thermal 
mass of Y induced by the coupling AY<f><l> c , which is of 
0(|A| 2 T 2 ) for |AY| < T. 

For simplicity, we will assume that all the dimension¬ 
less parameters appearing in the superpotential and the 
SUSY breaking scalar masses are of order unity. How¬ 
ever it should be noted that these parameters can have 
a variation of 0(0.1 10) easily. In particular, the super¬ 
potential parameters K n can have a much wider varia¬ 
tion without invoking fine-tuning. This gives us a rather 
large room to get an enough suppression of the axion an¬ 
gle fluctuation 56 through a stronger QCD before the /i- 
transition. At any rate, assuming that c X y > 0, t; X > 0 
and < 0 , the scalar potential (OBI) indeed yields the 
desired /x-transition as 

X = Y = 0 at t 

X ~ Y ~ y/ msuSY M P i at (28) 

with T(t,j) ~ tosusy- 

Now the Higgs and slepton fields can have a nontrivial 
evolution along the following flat direction: 

H d = (</>d, 0), L t = (<fig, 0), 

Hl = ( 0,V\faF+W), (29) 

which satisfies the F and D flat conditions. The relevant 
terms of the scalar potential of <f> d / are given by 

V 2 = TO?|<(>i | 2 + |</>i| 2 (HMi + h.c.) 

+ ( Eltf )(«£ +h ,.) 

+ (y~! 10j| 2 ) (l^rf | 4 + 4|0d ^| 2 + \(j)t \ A ), (30) 
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for fi = K\X 2 /Mpi, where 


m td — Cd,H 2 + gusY + 2|/i| 2 , 

m 4> e = c eH 2 + & to susy + M 2 - (31) 

Again, assuming c d ,e < 0 and &,£ < 0, but rn^ d ^ e (t 0 ) > 0 
due to //(to) ~ tosusy, the above scalar potential yields 

fa(tl) ~ ~ y/H(tl)Mpi, 

fa(t,j,) ~ 4>d,e(tfj,) ~ \/ msusyMpi, (32) 

and 

^d(to) = 0(100) GeV, &(i 0 ) = 0, 

fa (to) ~ A'(to) ~ Y(to) ~ V WSUSY Mpi. (33) 


To summarize, under a reasonably plausible assump¬ 
tion on the SUSY breaking during the primordial infla¬ 
tion and in the present Universe, the model with the su¬ 
perpotential (1251) can successfully realize the desired cos¬ 
mological evolution of the three relevant scales: the axion 
scale, the weak scale, and the QCD scale as given by (|28|) . 
© and (EH). Being generated by SUSY breakin g, an 
inflationary axion scale f a (ti) ~ ViL(f/)/msuSY fa (to) 
is determined to be well above the present axion scale 
f a (to) ~ V m SVSY Mpi, and a stronger QCD in the early 
Universe is realized to yield an enough suppression of 
the axion angle fluctuation even when H(tj ) saturates 
its upper bound. We note that the minimum of the ax¬ 
ion potential induced by the stronger QCD depends on 
arg(/t 3 A 3 ), but not on arg(R/j), while the minimum of 
the axion potential at present depends on arg(R//), but 
not on arg(K 3 A 3 ). As a result, the stronger QCD gen¬ 
erates an axion misalignment angle 0 m i s = 0(1), so that 
the axion dark matter with f l a = 12 dm arises naturally 
in our scheme. 

There is a remaining issue which should be addressed 
to complete our scheme. As we have noticed, the //- 
transition is foregone by a late-time thermal inflation. 
This suggests that the scheme should be accompanied by 
a late-time baryogenesis operating after the //-transition. 
In fact, the model of (l25l) offers an elegant mechanism 
to generate the baryon asymmetry through the rolling 
flat direction LH U |2l|. More detailed cosmology of our 
scheme, including the leptogenesis by rolling LH U , will 
be discussed elsewhere [221. 
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